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ABSTRACT 

With the release of the first year Fermi catalogue, the number of blazars detected above 100 
MeV lying at high redshift has been largely increased. There are 28 blazars at z > 2 in 
the "clean" sample. All of them are Flat Spectrum Radio Quasars (FSRQs). We study and 
model their overall spectral energy distribution in order to find the physical parameters of 
the jet emitting region, and for all of them we estimate their black hole masses and accretion 
rates. We then compare the jet with the accretion disk properties, setting these sources in the 
broader context of all the other bright 7-ray or hard X-ray blazars. We confirm that the jet 
power correlates with the accretion luminosity. We find that the high energy emission peak 
shifts to smaller frequencies as the observed luminosity increases, according to the blazar 
sequence, making the hard X-ray band the most suitable for searching the most luminous and 
distant blazars. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 



1 INTRODUCTION 

The Large Area Telescope (LAT) onboard the Fermi satellite de- 
tected, after 11 months of all sky survey, more than 1,400 sources, 
presented in Abdo et al. (2010a), with roughly half of them being 
BL Lac objects or Flat Spectrum Radio Quasars (FSRQs) (Abdo 
et al. 2010b, hereafter A10) and a few radio-galaxies, starbursts, 
galaxies, and Narrow Line Seyfert 1 galaxies. The corresponding 
catalog of AGNs at high Galactic latitude (|6| > 10°) is called 
First LAT AGN Catalog (1LAC). BL Lacs and FSRQs are approx- 
imately present in equal number. 

With respect to the previous sample (LAT Bright AGN Sam- 
ple, hereafter LBAS), constructed after 3 months of survey (Abdo 
et al. 2009; hereafter A09), the number of detected blazars is about 
7 times larger, as a result of the lower limiting sensitivity, ob- 
tained with the longer exposure and the smaller required signifi- 
cance (from 10 a of the first 3 months to the current 4 a level). 
Correspondingly, also the number of high redshift blazars detected 
in 7-ray s increased: in the LBAS there were 5 blazars at z > 2 
(and none at z > 3), while in the 1LAC catalogue there are 28 
sources at z > 2 (and 2 at z > 3). 

The increased number of high redshift 7-ray blazars allows 
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us to characterize them in a meaningful way, through their Spectral 
Energy Distributions (SEDs) and their modelling. Indeed, the cov- 
erage at other frequencies (besides the FermifLAT band) includes 
observations by the Swift satellite for all sources both in the optical- 
UV band (through the Optical-UV Telescope UVOT) and the soft 
X-ray band (0.3-10 keV, through the X-Ray Telescope XRT). 

It is also interesting to compare the properties of the high red- 
shift blazars detected in 7-rays with the high-z blazars detected 
in hard X-rays by the Burst Alert Telescope (BAT) instrument on- 
board the Swift satellite. All blazars at z > 2 are FSRQs, so, up 
to now, high redshift "blazars" coincide with high redshift FSRQs, 
since no BL Lac objects with a measured redshifts z > 2 has been 
detected so far. There are 10 FSRQs at 2 > 2 and 5 at z > 3 in 
the 3-year BAT all sky survey presented by Ajello et al. (2009), 
that have been studied in Ghisellini et al. (2010a, hereafter G10). 
The BAT and the LAT samples of high redshift blazars are rather 
well defined, since the sky coverage is quasi-uniform (excluding 
the Galactic plane) and we can consider these samples as flux lim- 
ited. 

The main aims of the present paper are then to characterize 
the properties of blazars detected at high energies at redshift greater 
than 2 and to see if we can understand the differences (if any) be- 
tween the blazars detected in the two bands (7-rays and hard X- 
rays). In G10, in fact, we suggested that the best way to select the 
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most powerful blazars at large redshifts is through a survey in the 
hard X-ray band, rather than in the 7-ray one, but this was based 
on small numbers. None of the 10 BAT blazars at z > 2 is present 
in the LB AS catalogue, and only 4 of them are in the 1LAC sample, 
i.e. have been detected after 1 1 months of survey by the FermifLKT 
instrument. In G10 we explained this through a change of the av- 
erage SED when the bolometric luminosity changes: by increasing 
it, the high energy hump of the SED peaks at smaller frequencies, 
in agreement with the blazar sequence as put forward by Fossati et 
al. (1998) and Donato et al. (2001), and interpreted in Ghisellini et 
al. (1998). 

We anticipate that our earlier suggestion remains valid, with 
important implications on the planned future hard X-ray survey 
missions, such as EXIST. 

In this paper we use a cosmology with h = Q,a = 0.7 and 
Qm = 0.3, and use the notation Q = W x Q x in cgs units (except 
for the black hole masses, measured in solar mass units). 



2 THE HIGH REDSHIFT SAMPLE 

We consider all blazars detected during the first year all-sky survey 
of Fermi and classified as "clean" in the catalogue of A10. These 
are all the blazars with |6| > 10°, detected at more than the 4a 
level whose identification is secure and unique. In total the 1LAC 
clean sample contains 599 sources (A10), of which 248 are FSRQs, 
all with a measured redshift, and 275 BL Lacs (116 with the red- 
shift measured). Among these, we selected the 27 blazars at z > 2 
as listed and classified by A10, plus an additional source, SWIFT 
J1656.3-3302 (z = 2.4), that Ghirlanda et al. (2010) recently clas- 
sified as FSRQ among the unidentified 1LAC sources. 

Five of these were already present in the LBAS list, i.e. the 
blazars detected at more than the 10cr level during the first 3 months 
of Fermi survey (A09). Four additional sources are present in the 
3-years survey of SwiftfBAT (A09) and they too have been studied 
in G10. 

Table Q] lists all sources: the top 19 blazars are studied in this 
paper, while the bottom 9 are the sources already present either in 
the BAT or LBAS samples. In this paper we present the spectral 
energy distributions (SED) and the modelling for the "new" ones, 
i.e. blazars not present in our previous study (G10). 



3 SWIFT OBSERVATIONS AND ANALYSIS 

For all blazars studied in this paper there are Swift observations. 
Even when they were performed during the 1 1 months of the 1LAC 
survey, they correspond to a "snapshot" of the optical-X-ray state 
of the source, while the 7-ray data are an average over the 11 
months. Given the very rapid blazar variability, the SEDs con- 
structed in this way should be considered, in all cases, not simulta- 
neous (but the Swift UVOT and XRT data are indeed simultaneous). 

The data were screened, cleaned and analysed with the soft- 
ware package HEASOFT v. 6.8, with the calibration database up- 
dated to 30 December 2009. The XRT data were processed with the 
standard procedures (XRTP IPELINE v. . 12 . 4). All sources 
were observed in photon counting (PC) mode and grade 0-12 (sin- 
gle to quadruple pixel) were selected. The channels with energies 
below 0.2 keV and above 10 keV were excluded from the fit and 
the spectra were rebinned in energy so to have at least 20-30 counts 
per bin in order to apply the x 2 test. When there are no sufficient 
counts, then we applied the likelihood statistic in the form reported 
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Table 1. List of blazars at 2 > 2. The upper part of the table reports the 
blazars studied in this paper (denoted by "0" in the last column); while the 
bottom part reports blazars studied previously: in Ghisellini et al. (2010; 
G10) (blazars with z > 2 detected by Swift/BAT); and in Ghisellini, 
Tavecchio & Ghirlanda (2009; G09) (blazars with z > 2 in LBAS with 
L 7 > lO^ergs" 1 ). 

by Cash (1979). Each spectrum was analysed through XSPEC v. 
12.5. In with an absorbed power law model with a fixed Galactic 
column density as measured by Kalberla et al. (2005). The com- 
puted errors represent the 90% confidence interval on the spectral 
parameters. Tab. [2] reports the log of the observations and the best 
fit results of the X-ray data with a simple power law model. The 
X-ray spectra displayed in the SED have been properly rebinned to 
ensure the best visualization. 

UVOT (Roming et al. 2005) source counts were extracted 
from a circular region 5"— sized centred on the source position, 
while the background was extracted from a larger circular nearby 
source-free region. Data were integrated with the uvotimsum 
task and then analysed by using the uvot source task. The ob- 
served magnitudes have been dereddened according to the formulae 
by Cardelli et al. (1989) and converted into fluxes by using standard 
formulae and zero points from Poole et al. (2008). Tab. [3] lists the 
observed magnitudes in the 6 filters of UVOT. 



4 MODELLING THE SED 

To model the SEDs of the blazars in this sample we used the same 
model used in G10. It is a one-zone, leptonic model, fully dis- 
cussed in Ghisellini & Tavecchio (2009). In that paper we em- 
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Table 2. Summary of XRT observations. The observation date column indicates the date of a single snapshot or the years 
during which multiple snapshots were performed. The corresponding note reports the complete set of observations integrated. 
The column "Exp" indicates the effective exposure in ks, while Nn is the Galactic absorption column in units of [10 20 cm -2 ] 
from Kalberla et al. (2005). T is the photon index of the power law model [F(E) oc E~ r ], Fq%— lOkoV ' s 0Dserve d 
(absorbed) flux. The two last columns indicate the results of the statistical analysis: the last column contains the degrees of 
freedom, while the last but one column displays the reduced x 2 or me value of the likelihood (Cash 1979), in the case there 
were no sufficient counts to apply the \ 2 test. 

a sum of observations of: 02/07/2007, 10/01/2008, 16/02/2008, 16/08/2009. 

b sum of observations of: 22/03/2007, 10/04/2007, 08/08/2007, 10/08/2007, 13/01/2008, 20/04/2008, 22/04/2008, 
26/10/2008, 22/04/2009. 

c sum of observations of: 26/12/2008, 08/01/2009. 
d sum of observations of: 21/02/2010, 25/02/2010. 
e sum of observations of: 06/03/2009, 01/10/2009, 04/10/2009. 

J Flux derived by using WebPIMMS with a rate of (5 ± 2) X 10 -3 c/s and fixed parameters. 
9 sum of observations of: 09/06/2006, 13/06/2006. 
h sum of observations of: 08/01/2009, 13/01/2009. 
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Table 3. Summary of Swift/UVOT observed magnitudes. Lower limits are at 3cr level. 
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Figure 1. SEDs of 0106+01 (=4C+01.02) together with the fitting mod- 
els, with parameters listed in Tab. [4] De-absorbed UVOT, XRT and BAT 
data are indicated by darker symbols (red in the electronic version), while 
archival data (from NED) are in light grey. Diamonds (and lower arrows, 
cyan in the electronic version) indicate UVOT data not de-absorbed by in- 
tervening Lyman-a clouds. The short-dashed line is the emission from the 
IR torus, the accretion disk and its X-ray corona. The solid thin (green) line 
is the non-thermal emission (sum of synchrotron, SSC and EC). The long 
dashed and the dot-dashed grey lines are the synchrotron self-Compton 
(SSC) and the External Compton (EC) components, respectively. The thick 
solid (blue) line is the sum of all components. The grey stripe in the 7-ray 
band corresponds to the FermifLfH sensitivity of the first 3 months (10cr, 
upper boundary) and for 1 1 months (4<r, lower boundary). 



phasize the relative importance of the different sources of the seed 
photons for the inverse Compton scattering process, and how they 
change as a function of the distance of the emitting region from the 
black hole. Here we briefly summarize the main characteristics of 
the model. 

The source is assumed spherical (radius R) and located at a 
distance i?di ss from the central black hole. The emitting electrons 
are injected at a rate Q(y) [cm -3 s _1 ] for a finite time equal to the 
light crossing time R/c. The shape of Q(y) we adopt is assumed 
to be a smoothly broken power law with a break at 7b: 



0(7) 
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1 + (7/7b)- sl+a 



(1) 



The emitting region is moving with a velocity /3c correspond- 
ing to a bulk Lorentz factor V. We observe the source at the viewing 
angle 9 V and the Doppler factor is 5 = l/[r(l — /3cos# v )]. The 
magnetic field B is tangled and uniform throughout the emitting 
region. We take into account several sources of radiation externally 
to the jet: i) the broad line photons, assumed to re-emit 10% of 
the accretion luminosity from a shell-like distribution of clouds lo- 
cated at a distance Rblb. = 10 17 L^^ 5 cm; ii) the IR emission 

from a dusty torus, located at a distance Rm. = 2.5 x W ls L^/^ 5 
cm; iii) the direct emission from the accretion disk, including its 
X-ray corona. Furthermore we take into account the starlight con- 
tribution from the inner region of the host galaxy and the cosmic 
background radiation, but these photon sources are unimportant in 
our case. All these contributions are evaluated in the blob comov- 
ing frame, where we calculate the corresponding inverse Compton 
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Figure 2. SED of PMN 0157-4614, B2 0242+23 and TXS 0322+222. Sym- 
bols and lines as in Fig. [2] 



radiation from all these components, and then transform into the 
observer frame. 

We calculate the energy distribution AT (7) [cm ] of the emit- 
ting particles at the particular time R/c, when the injection process 
ends. Our numerical code solves the continuity equation which in- 
cludes injection, radiative cooling and e pair production and re- 
processing. Ours is not a time dependent code: we give a "snap- 
shot" of the predicted SED at the time R/c, when the particle distri- 
bution N(y) and consequently the produced flux are at their maxi- 
mum. 

For all sources in our sample, the radiative cooling time of the 
particles is short, shorter than R/c even for low energetic particles. 
In Tab. [4] (last column) we have listed the values of j c , that is the 
minimum value of the random Lorentz factor of electrons cooling 
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Figure 3. SED of PKS 0420+022, PKS 0451-28, PKS 0458-02 and PKS 
0601-70. Symbols and lines as in Fig. [2] 
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SED of PMN 0625-5438, TXS 0907+230, TXS 0908+416 
-084. Symbols and lines as in Fig. [2] 



and 
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Figure 5. SED of TXS 1343+451, PMN 1344-1723, [WB92] 1537+2754 
and SWIFT J1656-3302. Symbols and lines as in Fig. El 
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Figure 6. SED of PMN 1959^1246, TXS 2118+188 and PMN 2135- 
Symbols and lines as in Fig. El 
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in one light crossing time. Since it is always smaller than 7b, al- 
most all the energy injected in the form of relativistic electrons is 
radiated away. Most of the cooling is due to the inverse Compton 
scattering with broad line photons, with a minor contribution from 
the synchrotron and the self-Compton process. Therefore we al- 
ways are in the fast cooling regime (i.e. 7 C < 7b). In this regime 
the produced luminosity does not depend on the amount of the radi- 
ation energy density, but only on the energy content of the injected 
relativistic electrons. 

Another implication is that, at lower energies, the AT (7) distri- 
bution is proportional to 7~ 2 , while, above 7b, N(y) oc 7~( S2+1 '. 
The electrons emitting most of the observed radiation have energies 
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7 P oak which is close to 7b (but these two energies are not exactly 
equal, due to the curved injected spectrum). 

The accretion disk component is calculated assuming a stan- 
dard optically thick geometrically thin Shakura & Sunjaev (1973) 
disk. The emission is locally a black body. The temperature profile 
of the disk is given e.g. in Frank, King & Raine (2002). Since the 
optical-UV emission is the sum of the accretion disk and the jet 
non-thermal components, for a few sources there is some degener- 
acy when deriving the black hole mass and the accretion rate. 

We model at the same time the thermal disk (and IR torus) ra- 
diation and the non-thermal jet-emission. The link between these 
two components is given by the amount of radiation energy den- 
sity (as seen in the comoving frame of the emitting blob) coming 
directly from the accretion disk or reprocessed by the BLR and the 
IR torus. This radiation energy density depends mainly on Pdiss, 
but not on the adopted accretion rate or black hole mass (they are 
in any case chosen to reproduce the observed thermal disk lumi- 
nosity). 

By estimating the physical parameters of the source we can 
calculate the power that the jet carries in the form of radiation (P r ), 
magnetic field (Pb), relativistic electrons (P c ) and cold protons 
(P p ) assuming one proton per electron. These powers are calcu- 
lated according to: 



P 



ttR 2 F 2 cU[ 



(2) 



where U' is the energy density of the i t h component in the comov- 
ing frame. 



4.1 Intervening Lyman-a absorption 

Being at z > 2 the optical-UV flux of the blazars in our sample 
could be affected by absorption of neutral hydrogen in intervening 
Lyman-a absorption systems. To correct for this, we use the attenu- 
ation calculated in G10 specifically for the UVOT filters, illustrated 
in Fig. 3 of that paper. 

Full details of our calculation will be described in Haardt et 
al. (in preparation), together with a more refined treatment of the 
mean attenuation and its variance around the mean. The current 
procedure is very crude, especially when the attenuation is large 
(i.e. optical depths larger than unity) because in such cases most of 
the attenuation is due to very few clouds, implying a large variance. 
However, we note that the variance of the attenuation is largely 
reduced when the actual filter width is taken into account (Madau 
1995). Our absorption model results in a mean number of thick 
systems which is < 1 for z<4, so we do not expect excessive off- 
set of the attenuation along individual line of sight with respect to 
the mean value. 

When presenting the SED of our sources, we will show both 
the fluxes and upper limits de-reddened for the extinction due 
to our Galaxy and the fluxes (and upper limits) obtained by de- 
absorbing them with the r e g shown in Fig. 3 in G 10. 



5 RESULTS 

Table|4]lists all parameters used to model the SEDs of our blazars, 
Table [5] lists the different forms of power carried by the jet and 
Fig. HHU show the SEDs of the 19 blazars studied in this paper 
and the corresponding fitting model. In all figures we have marked 
with a grey shaded area the FennifLAT sensitivity, bounded on the 
bottom by considering one year of operation and a 5<r detection 
level, and on the top by considering 3 months and a 10<j detection 
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Figure 7. Distribution of the black hole masses derived for the z > 2 
FermrTLAT (top) and SwiftfBKt samples (bottom). The hatched area in the 
top panel corresponds to the 4 blazars in common. 



level (this assumes a common energy spectral index of a 7 ~ 1, the 
sensitivity limit for other spectral indices is slightly different, see 
Fig. 9 in A10). All these sources were not detected in the first 3 
months, in fact the (11 months) 7-ray data points are very close to 
the lower boundary of the grey area. There are exceptions: 0106+01 
(=4C+01.02) is brighter than the 3-months, 10ct sensitivity limits 
even if it has not been included in LB AS. This is due to a rather 
strong variability of the source, fainter in the first 3 months and 
brighter soon after. The same occurred for 1344—1723 and 0451- 
28. The opposite happened for 0227-369, 0347-21 1 and 0528+134 
(i.e. they were brighter during the first 3 months), but their flux, 
averaged over 1 1 months, was in any case large enough to let their 
inclusion in the 1LAC sample. 

Some of the sources have a sufficiently good IR-optical-UV 
coverage to allow to see a peak of the SED in this band (see 
for instance 0420+022; 0451-28; 0458-02; 0625-5428; 0907+230; 
0908+416; 1149-084; 1656-3302). The other sources have a SED 
consistent with a peak in this band, but the lack of data also al- 
lows for a peak at lower frequencies. We interpret the peak in the 
optical band as due to the accretion disk, and assume its presence 
also in those blazars where it is allowed, but not strictly required. 
By assuming a standard Shakura-Sunyaev (1973) disk we are then 
able to estimate both the black hole mass and the accretion rate. 
This important point has been discussed in G10, in Ghisellini et al. 
(2009) (for S5 0014+813) and in Ghisellini & Tavecchio (2009). 

The radio data cannot be fitted by a simple one-zone model 
specialized to fit the bulk of the emission, since the latter must be 
emitted in a compact region, whose radio flux is self-absorbed up 
to hundreds of GHz. The radio emission should come from larger 
regions of the jet. On the other hand, when possible, we try to have 
some "continuity" between the non-thermal model continuum and 
the radio fluxes (i.e. the model, in its low frequency part, should not 
lie at too low or too high fluxes with respect to the radio data). In 
the following we briefly comment on the obtained parameters. 

Dissipation region — The distance Pdiss at which most of the dis- 
sipation takes place is one of the key parameters for the shape 
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Table 4. List of parameters used to construct the theoretical SED. Not all of them are "input parameters" for the model, because -Rblr is uniquely determined 
from L d , and the cooling energy -y c is a derived parameter. Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10 cm and (in parenthesis) 
in units of Schwarzschild radii; Col. [4]: black hole mass in solar masses; Col. [5]: size of the BLR in units of 10 15 cm; Col. [6]: power injected in the blob 
calculated in the comoving frame, in units of 10 45 erg s — 1 ; Col. [7]: accretion disk luminosity in units of 10 45 erg s _1 and (in parenthesis) in units of Lecm; 
Col. [8]: magnetic field in Gauss; Col. [9]: bulk Lorentz factor at -Rdiss; Col. [10] and [11]: break and maximum random Lorentz factors of the injected 
electrons; Col. [12] and [13]: slopes of the injected electron distribution [Q("f)] below and above 7b; Col. [14] values of the minimum random Lorentz factor 
of those electrons cooling in one light crossing time. The total X-ray corona luminosity is assumed to be in the range 10-30 per cent of L d . Its spectral shape 
is assumed to be always oc v~ 1 exp(— hu/150 keV). The viewing angle 8 V is 3° for all sources. 



of the overall SED, since it controls the amount of energy den- 
sities as seen in the comoving frame (see Ghisellini & Tavecchio 
2009). For almost all sources we have i?di ss < ^blr, while for 
0106+01, 0907+230, 1343+451, 1344-1723, 1959^1246, R diss is 
slightly larger than i? B LR, and for 0743+259 R d i ss ~ 2i? B LR ■ 

In all sources the dominant cooling is trough inverse Comp- 
ton off the seed photons of the BLR. This is true also for the 
few blazars in which -Rdiss -Rblr since, even if (7 B lr seen 
in the comoving frame is smaller, also Ub is smaller, imply- 
ing that the ratio £/blr/^b is similar to the values in other 
sources [Ub LK /U b ranges between ~30 (1344-1723) and ~100 
(0106+01, 0907+230)]. However, the decreased cooling rate in 
0743+259 makes 7 C to be larger (7 C = 102, see Tab.©. 

With larger still -Rdiss 2> -Rblr, the main seed photons for the 
Compton scattering process would became the photons produced 
by the the IR torus (if it exists), but this case does not occur for our 
sources. 

Compton dominance — This is the ratio between the luminosity 
emitted at high frequencies and the synchrotron luminosity. The 
average magnetic field is found to be of the order of 1 Gauss, with 



a corresponding magnetic energy density that is around two orders 
of magnitude lower than the radiation energy density. Correspond- 
ingly, all sources are Compton dominated. 

Black hole masses — Fig. [7] shows the distribution of black hole 
masses for the 28 blazars at z > 2 and compares them with the dis- 
tribution of masses for the high redshift BAT blazars. Although the 
black hole masses of the BAT sample extend to larger values, there 
are still too few sources to estimate if the two distributions are dif- 
ferent. It is interesting to note that all but 3 sources (0420+022, 
0907+230 and 2135-5006) have black hole masses greater than 
1O 9 M . In Ghisellini, Tavecchio & Ghirlanda (2009) we consid- 
ered the Fermi blazars of 7-ray luminosity L 7 > 10 48 erg s _1 , 
finding, with the same method and model applied here, that for all 
these sources the black hole mass was greater than a billion solar 
masses. Therefore all blazars with L 7 > 10 48 erg s _1 have black 
holes heavier than 1O 9 M0, while the vast majority, but not all, 
blazars at z > 2 have such large black hole masses. We searched 
in the literature other estimates of the black hole masses for the 
objects in this sample, finding M = 2.3 x 1O 9 M for 0836+710 
(estimated by Liu, Jiang & Gu 2006), and other few limits for the 
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Table 5. Logarithm of the jet power in the form of radiation (P r ), Poynting 
flux (Pb), bulk motion of electrons (P c ) and protons (P p , assuming one 
proton per emitting electron). Powers are in erg s — . The bottom part of the 
table reports the data derived in G10 and G09. 



black hole masses for 0836+710 and 0528+134, that were however 
based assuming an isotropic 7-ray emission. 

Disk luminosities — We are considering very powerful blazars, so 
we do expect large disk luminosities, not only on the basis of an 
expected positive trend between the observed non-thermal (albeit 
beamed) and the accretion luminosities, but also on the basis of 
the observed luminosities of the broad lines, that should linearly 
depend on the accretion power. What is interesting is that all the 
FSRQs analyzed up to now (i.e. belonging to the LBAS sample or 
to the subset of high redshift 1LAC and BAT samples) have a ratio 
-kd/^Edd between 10~ 2 and 1. This can be seen in the mid panel of 
Fig. [8] that shows La as a function of the derived black hole mass. 
The two dashed lines correspond to the Eddington and 1 % Edding- 
ton luminosities. This confirms the idea of the "blazars' divide" as a 
result of the changing of the accretion mode (Ghisellini, Maraschi 
& Tavecchio 2009): from the standard Shakura-Sunyaev (appro- 
priate for all FSRQs) to the ADAF-like regime (appropriate for BL 
Lacs). The z > 2 blazars analysed here have I/d/^Edd ratios rang- 
ing from 0.05 and 0.7. The exact values of the disk luminosities 
derived here are the frequency integrated bolometric luminosities 
of the assumed Shakura-Sunyaev accretion disk model that best 
interpolates the data. On the other hand, any other accretion disk 
model has to fit the data as well, implying that our values of La 
are robust, and nearly model-independent, within the limit of the 
uncertainties of the observed data. 
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Figure 8. The observed 7-ray luminosity in the 0.1-10 GeV band, the 
accretion luminosity and the total jet power Pj ct as a function of the de- 
rived black hole mass. All points correspond to FSRQs. Different symbols 
correspond to LBAS FSRQs with L~ f > 10 48 erg s _1 (red stars, analysed 
in Ghisellini, Tavecchio & Ghirlanda 2009); 1LAC FSRQs with z > 2 
(empty blue circles); BAT FSRQs with 2 > 2 (black diamonds, G10) and 
the LBAS FSRQs with L 7 < 10 48 erg s — 1 (grey squares, Ghisellini et al. 
2010b). The 7-ray luminosity (top panel) is the observed beamed one, and 
it has not to be confused with P r , the power spent by the jet to produce the 
radiation we see. The mid panel shows that all FSRQs have disk luminosi- 
ties between 0.01 and 1 Eddington luminosity. The bottom panel shows that 
Pjet can be larger (but not by a big factor) than the Eddington luminosity 
corresponding to the dashed line. 

Jet powers — The values listed in Tab.|5]are very similar to the val- 
ues derived for other powerful Fermi FSRQs. They are not, how- 
ever, the absolutely greatest powers found. This can be seen in the 
bottom panel of Fig. [8] showing Pj c t as a function of the black hole 
mass, and in the mid and bottom panels of Fig. [9] where we plot the 
power of the jet spent in the form of radiation (P r ) and the total jet 
power Pj c t as a function of id. We can compare the z > 2 1LAC 
FSRQs with those present in the LBAS catalogue and the high red- 
shift BAT blazars. Not surprisingly, we see that in these planes the 
1LAC z > 2 sources follow the distribution of the most luminous 
7-ray blazars. Remarkably though, the z > 2 BAT FSRQs appear 
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Figure 9. Top: the observed 7-ray luminosity L 1 as a function of the 
accretion luminosity Ld for the LBAS FSRQs (grey squares), the 2 > 2 
1LAC (blue circles) and Svvi/f/BAT (diamonds) FSRQs. The dashed line 



indicates L~. 



L,j. Mid: the power P r as a function of L^. P r can be 



considered as a robust and almost model-independent lower limit to the jet 
power. Bottom: the total jet power Pj ct = Pb + P c + Pp as a function 
of L,j. Almost all sources have Pj ct > L^. One cold proton per emitting 
electron is assumed. 



to lie at the extreme of the distributions, being the more powerful 
in Ld, and among the most powerful in P r and Pj ct . For calculat- 
ing the power carried by the jet in the form of protons, we re-iterate 
that we have assumed one cold proton per emitting electron: if there 
exist a population of cold electrons, and no electron-positron pairs, 
than we underestimate P p and then Pjct, while if there are no cold 
leptons but there are pairs then we overestimate P p . Finally, pro- 
tons are assumed cold for simplicity (and "economy"), but they 
could be hot or even relativistic (if, e.g. shocks accelerate not only 
electrons but also protons), and in such cases the power is underes- 
timated. For a detailed discussion about the presence of electron- 
positron pairs in blazars' jets we refer to the discussions in Sikora 
& Madejski (2000), Celotti & Ghisellini (2008) and G10, where 
one can finds arguments limiting the amount of pairs in the jet. A 
few electrons-positrons per proton are possible, but not more. 

Jet powers vs accretion luminosities — Fig.|9]shows that the corre- 
lations found in G10 between P r and/or Pj ct and Ld are confirmed. 
We remind the reader that P r and Ld are independent quantities 
even if the main radiation mechanism is the inverse Compton pro- 
cess using broad line photons as seeds, that in turn are proportional 
to the accretion disk luminosity. This is because the radiative cool- 
ing of the emitting electrons is complete, implying that the pro- 
duced jet luminosity becomes independent on the amount of radi- 
ation energy density. In other words: in the fast cooling regime the 
jet always emits all the energy of its relativistic electrons, no matter 
the amount of the luminosity of the accretion disk. 

A least square fit returns a chance probability P = 4 x 10~ 8 
that log P r and log Pj ot are correlated with log Ld (and the cor- 
relation are consistent with being linear). They remain significant 
also when considering the common redshift dependence, although 
the chance probability increases to P = 4 x 10~ 4 (for the P r -Ld 
correlation) and to P = 10~ 3 (Pj ct -L d ). 

As expected, the 1LAC blazars at high redshifts are among 
the most powerful, even if there are blazars at lower redshifts with 
comparable powers. This can be seen comparing the empty circles, 
corresponding to the 1LAC blazars of our sample, with the LBAS 
FSRQs of L 7 > 10 48 erg s _1 (stars) and the BAT FSRQs at 2 > 2 
(black diamonds). There are a few sources with P r > Ld, and 
several with P r ~ Ld. The jet in these blazars, only to produce the 
radiation we see, requires a power comparable to (or even larger 
than) the disk luminosity. The P r power should be considered a 
very robust estimate of the minimum jet power: it is robust because 
it is almost model-independent (P r ~ L 7 /r 2 , see G10), and it is 
a lower limit because it corresponds to the entire jet power being 
converted into radiation at the 7-ray emitting zone. Indeed, if there 
is one proton per emitting electron, the total jet power, dominated 
by the bulk motion of cold protons, becomes a factor ~ 10 larger 
than Ld (bottom panel of Fig. |5J, with the FSRQs of our sample 
distributed in a large portion of the Pjct-Ld plane. 

We believe that the relation between both P r and Pj ct with Ld 
is a key ingredient to understand the birth of jets: accretion must 
play a key role. 

Comparison with other models — Several groups (Larionov et al. 
2008; Marscher et al. 2008, 2010; Sikora, Moderski & Madejski 
2008) proposed that the emitting region, especially during flares, is 
produced at distances from the central black hole of the order of 
10-20 pc (much larger than what we assume) at the expected loca- 
tion of a reconfinement shock (e.g. Sokolov, Marscher & McHardy 
2004). On the basis of an observed peculiar behaviour of the polar- 
ization angle in the optical, Marscher et al. (2008) thus suggested 
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Figure 10. Top: The ratio of the (rest frame) luminosities in the 0.1— 
100 GeV and 15-55 keV bands as a function of the rest frame peak fre- 
quency of the high energy bump. For sources not detected in one of the two 
bands, we have used the corresponding flux resulting from the modelling of 
their overall SED. Different symbols are from the sources detected only by 
Fermi/LAT, by both Swift/BAT and Fermi/LAT and only by Swift/BAT, as 
labelled. For comparison we also show all FSRQs (grey dots) and BL Lacs 
(squares) in the LBAS sample. The continuous line shows the estimate us- 
ing the smoothly broken power law function (see text) with a x = 0.5 and 
a 7 = 1.6. Bottom panel: the same, but plotting observed peak frequen- 
cies. The continuous lines are for (z) = 2.5 and (z) = 0.5. See how the 
simple function in Eq. [3] interpolates well the data of FSRQs. The X and 
7-ray SED of BL Lacs, instead, is not due to a single radiation process, 
since the hard X-rays are often due to the tail of the synchrotron emission. 
As a consequence, they show the opposite behavior of FSRQs: a smaller 
Lx/Lj ratio when increasing ^ pca k> indicating an increasingly stronger 
contribution of the synchrotron flux to the hard X-ray band. 



that blobs ejected from the central region are forced by the mag- 
netic field to follow a helical path, accounting for the observed ro- 
tation of the polarization angle in the optical. Flares (at all wave- 
lengths) correspond to the passage of these blobs through a stand- 
ing conical shock, triggered by the compression of the plasma in 
the shock. This has important consequences for the variability of 
the emission: since the emission region is located at large distances 
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Figure 11. The rest frame peak frequency f poa k of the high energy emis- 
sion as a function of the kpeak^poak luminosity. Symbols as in Fig. 1101 
A clear trend can be seen when considering BL Lacs and FSRQs together, 
while FSRQs only are characterized by a very large dispersion. High red- 
shifts FSRQs occupy the region of the largest luminosities and smallest 
^peak- We indicate the two "outliers": BL Lac and Ap Lib. 



from the central engine, its size is large, and the expected variabil- 
ity timescale cannot be very short. Assuming J?diss = 15 pc, a 
jet aperture angle of 6j ct 
variability time scale of t 
of 1.5 (1+z) months, that for sources at z > 2 implies a mini- 
mum variability timescales of 5-6 months. The main high energy 
emission mechanism is still the inverse Compton process, using 
as seeds the IR radiation of a surrounding torus (Sikora, Moderski 
& Madejski 2008) with a possible important contribution from jet 
synchrotron radiation (Marscher et al. 2008, 2010). The main dif- 
ficulty of these models concerns the expected variability, predicted 
to occur on a very long time scales, if the size of the emitting re- 
gion is proportional (through e.g. the opening angle of the jet), to 
the distance of the source to the black hole. Instead the observed 
7-ray flux in all strong 7-ray sources (the ones for which a reliable 
variability behaviour can be established) varies on much shorter 
time scales, and factor 2 flux changes can occur even on 3-6 hours 
(see Tavecchio et al. 2010 for 3C 454.3 and PKS 1510-089; Bon- 
noli et al. 2010; Foschini et al. 2010 and Ackermann et al. 2010 
for 3C 454.3; Abdo et al. 2009b for PKS 1454-354; Abdo et al. 
2010c for PKS 1502+105). This indicates that the source is com- 
pact. This in turn suggests (although it does not prove) that its lo- 
cation cannot be too far from the black hole. It then also suggests 
that it is within the broad line region. In turn, this suggests that 
the broad lines are the main seeds for the inverse Compton scatter- 
ing process. Occasionally, though, dissipation could occur further 
out, where the main seeds are the infrared photons produced by a 
putative torus surrounding the accretion disks. Since the seed pho- 
tons have smaller frequencies, in these cases the produced high en- 
ergy spectrum suffers less from possible effects of the decreasing 
(with seed frequencies) scattering Klein-Nishina cross section and 
less from possible photon-photon interactions leading to electron- 
positron pair production. The decreased importance of both effects 
may account for high energy spectra extending, unbroken, up to 
hundreds of GeV. These cases should be characterized by a longer 
variability timescale. 

In the model of Marscher et al. (2008) a very short variability 
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timescale still indicates a very compact emitting region, but never- 
theless located at a large distance from the central black hole. 



5.1 Comparing BAT and LAT high redshift blazars 

Both the 1LAC sample and the BAT 3-years survey have a rather 
uniform sky coverage, and both approximate a flux limited sample. 
The 1LAC sample has a limiting flux sensitivity that depends on the 
7-ray spectral index of the sources, but since we are dealing with 
FSRQs only (whose a 7 is contained in a relatively narrow range), 
we can consider this sample as flux limited. It is then interesting 
to compare the high redshift blazars (all of them are FSRQs) con- 
tained in the two samples. We remind that for blazars with z > 2, 
we have 10 FSRQs in the BAT sample, 28 in the 1LAC, and 4 in 
both. The X-ray to 7-ray SED of these sources is very similar: 
even if we have only 4 blazars in common, for all sources the SED 
has a high energy peak in the ~MeV-100 MeV band, with a x < 1 
and a 7 > 1. 

Thus the spectrum can be approximated by a broken power 
law. For illustration, consider the smoothly broken power law of 
the form 



L(u) tx 



(f/"b)~ a " 



1 + (v/vb) a ~r- 



(3) 



If the energy indices a x < 1 and a 7 > 1, the peak is at i/ pca k = 
^b[(l - a x )/(pLy - l)] 1 /^!"^). With this function we can easily 
calculate the ratio of the BAT [15-55 keV] to LAT [0.1-100 GeV] 
luminosities as a function of i/ pea k, and see if it compares well to 
the data. We alert the reader that by "data" we mean real observed 
data when the source has been detected in the X-ray and 7-ray 
band, while, when the detection is missing, we mean the "data" 
coming from our fitting model described in §4 (and in more detailed 
in Ghisellini & Tavecchio 2009). 

The ratio Lx / L- t as a function of f pea k calculated using Eq. 
[3]setting a x — 0.5 and a 7 = 1.6 is plotted in Fig.[l0]as a grey (or- 
ange in the electronic version) line, together with the points corre- 
sponding to high redshift blazars, studied in this paper and in G10. 
Furthermore, Fig.[l0]reports also the data of all the LBAS blazars 
studied in Ghisellini et al. (2010b). These are FSRQs (filled grey 
circles) and BL Lacs (empty grey squares). In the top panel ^ pca k 
is in the rest frame of the source, while in the bottom panel ^ pca k 
is the observed one. 

Fig.[l0]shows a remarkable agreement between the FSRQ data 
and the simple broken power law of Eq.[3] both for high redshifts 
and for less distant FSRQs. BL Lac objects, instead, are not well 
represented by Eq.[3] In fact, in many BL Lacs the hard X-rays cor- 
respond to the (steep) tail of the synchrotron component. Therefore 
the X and the 7-rays are produced by two different mechanisms, 
and Eq.[3]does not represent their overall high energy SED. In BL 
Lacs the importance of X-rays increases increasing ^ pca k, because 
of the increasing importance of synchrotron flux in the hard X-rays. 

Coming back to high redshift FSRQs, only 4 sources have de- 
tection in both bands, while 24 FSRQs are detected only by the 
LAT and 6 only by the BAT instrument. Our model explains the 
large fraction of sources that are detected only in one of the two 
instruments as due to the different ^ pca k of the sources: if ^ pca k 
is large (above 10 MeV), Fj/Fx is large and the source is rela- 
tively weak in the BAT band, while if f pe ak <10 MeV the source 
becomes relatively weak in the LAT band and strong in the BAT. 

It is interesting to see if the derived f pe ak correlates with 
the bolometric luminosity. Fig. [TT] shows that a trend indeed ex- 
ists: more powerful sources have smaller i-'peak (we have used 



VpeakLv eak as a proxy for the bolometric luminosity). But com- 
paring with all the LBAS FSRQs (grey dots) we see that the z > 2 
luminous FSRQs belong to a broader (i.e. more scattered) distribu- 
tion, and that the high-z BAT blazars are really the most extreme. 

We can conclude that i) there is a a trend between the high 
energy peak and the peak luminosity, ii) that this correlation has a 
large scatter, even if iii) the z > 2 FSRQs show the same trend 
with less scatter (but this may be due to the still small number) and 
finally iv) the z > 2 blazars in the 3-years BAT survey all lie in the 
highest luminosity, smallest u pca ^ part of the plane. 

When more BAT detections of high redshift LAT blazars will 
become available (and, conversely, when LAT will detect more 
BAT high-z blazars) this trend can be tested directly (i.e. with- 
out modelling the SED). The importance of this is two-fold: first 
we can estimate in a reasonable way the peak energy of the high 
energy emission having the hard X-ray and the 7-ray luminosities, 
and second (and more important) we could conclude that the most 
powerful blazars can be more easily picked up trough hard X-ray 
surveys, as the one foreseen with the EXIST mission (Grindlay et 
al. 2010). 



6 SUMMARY AND DISCUSSION 

The total number of z > 2 blazars with high energy information 
is 34 (28 with a LAT detection, 6 with BAT, and 4 with both). It is 
still a limited number, the tip of the iceberg of a much larger (and 
fainter) population, but it is derived from two well defined samples 
(LAT and BAT), that we can consider as flux limited and coming 
from two all sky surveys (excluding the galactic plane). The main 
result of studying them is that all the earlier findings concerning 
the physical parameters of the jet emitting zone, the jet power, and 
the correlation between the jet power and the disk luminosity are 
confirmed. 

They are in agreement with the blazar sequence, i.e. their non- 
thermal SED are "redder" than less luminous blazar, with a large 
dominance of their high energy emission over the synchrotron one. 
This implies that the disk emission is left unhidden by the syn- 
chrotron flux, and this allows an estimate of the black hole mass 
and the accretion rate. The uncertainties associated with these es- 
timates are relatively small within the assumption that the thermal 
component is produced by a standard Shakura-Sunyaev disk with 
an associated non-spinning hole. In G10 we argued that in any case 
the masses are not largely affected by this assumption, and in par- 
ticular that in the Kerr case the derived masses are not smaller, de- 
spite the greater accretion efficiency. The possibility of an intrin- 
sic collimation of the disk radiation appears more serious. If the 
disk is not geometrically thin, but e.g. a flared disk, then we ex- 
pect a disk emission pattern concentrated along the normal to the 
disk, i.e. along the jet axis. We argued previously (Ghisellini et al. 
2009) that this can be the case of S5 0014+813 at z = 3.366, 
an extraordinary luminous blazars (detected by BAT) with an esti- 
mated "outrageous" black hole mass of 4 x 10 10 M o . And indeed 
we found it to be an "outlier" in the P T -Ld plane. Reverting the 
argument, it implies that the other FSQRs, obeying a well defined 
P T -Ld trend, should have disk luminosities with a quasi-isotropic 
pattern, i.e. standard, not flaring, disks. The other severe uncertainty 
on the mass estimation for objects at large redshifts is the amount of 
attenuation of their optical-UV flux, due to intervening Lyman-a 
clouds. We have corrected for this assuming an average distribution 
of clouds, and when the attenuation is due to a few thick clouds 
the expected variance is large. While we plan to refine such esti- 
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Figure 12. The ratio between the luminosities in the BAT and LAT energy 
ranges (i.e. 15-55 keV and 0.1-10 GeV) as a function of the disk luminos- 
ity, in Eddington units. Big dots are the 6 BAT blazars in the z > 2 sample 
not detected (yet) by Fermi/LAT, squares are BAT blazars of the same sam- 
ple already present in the 1LAC sample, triangles are the sources discussed 
in this paper, and the small grey dots are all the FSRQs present in the LBAS 
sample discussed in G10. The dashed line is the best least square fit (chance 
probability P = 2 X 10" 6 ). 



(i.e. a blazar) there must exist ~ 2F similar sources pointing in 
other directions, the finding of even a few blazars at large red- 
shifts with a large black hole mass can put very interesting con- 
straints on the black hole mass function of jetted sources. This is- 
sue has been discussed by us previously (G10), and we suggested 
that the existence of the blazar sequence, plus an important K- 
correction effect, makes the hard X-ray range the best band where 
to search for the high-z heaviest black holes. Here we re-iterate 
this suggestion offering a supplementary information. For each an- 
alyzed FSRQs (belonging to the LBAS sample, or in this paper), 
we have calculated the ratio between the expected BAT luminos- 
ity [15-55 keV, rest frame] and the observed [0.1-100 GeV, rest 
frame] LAT luminosity. This ratio is not an observed quantity, since 
very few FSRQs (in comparison with LAT) have been detected by 
BAT, but it is a result of the model. Then, with the same model, 
we have calculated the disk luminosity in units of Eddington. Fig. 
1121 shows the Lbat/^lat ratio as a function of Ld/^Edd- We 
have marked with different symbols the FSRQs analyzed in this 
paper and the ones analysed previously (G10). We have excluded 
BL Lacs for which we have only an upper limit on the disk lu- 
minosity. Fig. [T2] suggests a clear trend (albeit with some scatter): 
FSRQs accreting close to Eddington emit relatively more in the 
hard X-rays than above 100 MeV. Formally, a least square fit re- 
turns (Z/bat/£bat) oc (Ld/^Edd) 1 ' 2 with a chance probability 
P = 2. x 10" 6 (for 95 objects). This implies that hard X-ray sur- 
veys benefit of a positive bias when looking for blazars with black 
holes accreting close to Eddington. In turn, at high-z, this implies 
the finding of the heaviest black holes, since it is very likely that at 
those early epochs (e.g. z > 2) all black holes are accreting close 
to the Eddington rate. 



mates (Haardt et al. in preparation), it is unlikely that this implies 
a systematic error on the mass estimates, leading to larger values: 
statistically, the mass distribution should not be seriously affected 
by this uncertainty. 

Although affected by the same issues outlined above, also the 
correlation between the jet power and the disk luminosity should 
resist when a more refined treatment of the Lyman-a absorption 
will be available, by the same arguments. Therefore we can con- 
clude that the accretion rate is really a fundamental player in pow- 
ering the relativistic jet, and we refer the reader to Ghisellini et al. 
(2010b) for a detailed discussion about this finding. If true, these 
ideas lead us to suggest that powerful, high redshifts "true" (i.e. 
really lineless) BL Lacs do not exist. 

Another, at first sight surprising, result of our study is that the 
correlation between the peak frequency of the high energy emis- 
sion and the 7-ray luminosity at this peak frequency exists also 
for the z > 2, highly luminous FSRQs. It is surprising because 
in very powerful sources the radiative cooling is complete (i.e. all 
electrons with 7 larger than an few cool in one light crossing time). 
Therefore the electrons responsible for the peak have energies that 
are not fixed by radiative cooling, but by the injection function [i.e. 
by the Q(j) function given in Eq. 1, and more precisely by the 
value of 7b]. On the other hand, when considering all blazars in 
the LBAS sample, we see that the scatter, for large luminosities, is 
much larger, so the apparent correlation for the z > 2 blazars can 
be due to small statistics. Clearly, it is a point to investigate further. 

It is interesting to ask what will be the best strategy to find the 
most luminous FSRQs at the largest redshifts. The interest lies in 
the link between jet power and disk luminosity: finding the most 
powerful jets implies to find the most accreting systems, hence 
the heaviest black holes. Since for each source pointing at Earth 



7 CONCLUSIONS 

We summarize here our main conclusions: 

• The blazars detected by Fermi at z > 2 are all FSRQs, with 
typical "red" SEDs. 

• These FSRQs are very luminous and powerful, but they are not 
at the very extreme of the distribution of luminosity and jet power. 

• These sources have heavy black holes (M ~ 1O 9 M0) and ac- 
cretion luminosities greater than ~10% Eddington. When includ- 
ing all FSRQs in the LBAS sample, irrespective of redshift, the 
accretion disk luminosities is greater than 1% Eddington. 

• The trend of redder SED when more luminous (i.e. one of 
the defining characteristics of the blazar sequence) is confirmed, 
and it is even present within the relatively small range of observed 
luminosity of the z > 2 blazars. 

• The correlation between the jet power and the disk luminos- 
ity is confirmed and points to a crucial role played by accretion in 
powering the jet. 

• FSRQs with accretion disks closer to the Eddington luminos- 
ity have jets emitting a "redder" SED, and therefore can be more 
efficiently picked up by hard X-ray surveys (such as the one fore- 
seen by EXIST), rather than by surveys in the hard 7-ray band. 
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